Chapter 3 

Linear Response Theory 



3.1 Measurements and Correlation Functions 

We will now discuss how to represent physical measurements. All local phys- 
ical measurements of a macroscopic quantum (many-body) system amount 
in practice to a process in which a localized disturbance is created by an 
applied external force, in the neighborhood of some point f at some time t, 
and the response of the system is then measured at some other point r ' at 
tome later time t' > t. On the other hand, in many cases we also will be 
interested in the global response of a system to an uniform perturbation, as 
in the case of most thermodynamic measurements. 

There are many examples of these type of measurements. For example, 
in a typical optical experiment one has an external electromagnetic wave im- 
pinging on a system (for instance, a metal). In the case of an electromagnetic 
field with sufficiently low energy we can ignore its quantum mechanical nature 
and treat it as a classical wave. It will interact with the degrees of freedom 
in a variety of ways. Thus, the scalar potential A couples to the local charge 
density through a term in the Hamiltonian of the form J al 3 xA (x,t)p(x,t), 
where p(x, t) is the second quantized charge density operator (in the Heisen- 
berg representation). On the other hand, the vector potential A(x, t) couples 
to the particle current operator j (x, t) (the full form of the coupling is actu- 
ally more complex, as we shall see below). Thus, at this level, we see that an 
external classical electromagnetic field can induce density waves and currents 
in a system. 

However, at higher energies the quantum nature of the electromagnetic 
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field cannot be ignored, and in reality what one is considering are the scatter- 
ing processes of external photons with the internal degrees of freedom of the 
system of interest. Photons couple to the atoms and through this coupling 
we can gain information on the phonon degrees of freedom. Thus, Raman 
scattering of photons provide a direct window to the dynamics of phonons 
(among many other things). 

On the other hand, X-ray photons are more energetic and sufficiently so 
that they can in fact eject an electron from the system, leaving behind a 
hole. Thus, by measuring the properties of the ejected photo- electron it is 
possible to infer the spectral properties of the hole. This is an example of 
a photoemission experiment. If the momentum and energy of the photo- 
electron are measured, this is an angle-resolved photo-emission (ARPES) 
experiment, which essentially measures the retarded Green function of the 
hole. 

Other scattering processes couple instead to the spin of the degrees of 
freedom. For example in neutron scattering experiments, the initial state has 
neutrons with well defined energy, momentum and spin polarization which 
scatter off the spin degrees of freedom into a final state also with well defined 
energy momentum and spin polarization. The leading coupling is due to 
the effectively short-range interaction between the spin of the neutron and 

— # 

the local spin density S(x, t) of the degrees of freedoms of the system. The 
interaction between the neutron magnetic moment and the local spin den- 
sity leads to a term on the Hamiltonian which is linear in S(x,t). Neutron 
scattering experiments measure the spin-spin correlation function. 

We will now discuss a general theory of these type of measurements. Let 
H be the full Hamiltonian describing the system in isolation. One way to 
test its properties is to couple the system to a weak external perturbation 
(assuming that the ground state is stable) and to determine how the ground 
state and the excited states are affected by the perturbation. Let 0(x, t) be 
a local observable, such as the local density, the charge current, or the local 
magnetization. The (total) Hamiltonian H T for the system weakly coupled 
to an external perturbation, which we will represent by a Hamiltonian H ext , 
is 

H T = H + H cxt (3.1) 

We will assume that the perturbation is not only weak, in which case 
its effects should be describable in perturbation theory, but it is also adia- 
batically switched on and off. Then as we discussed above, the Heisenberg 
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representation of the isolated system, with Hamiltonian H, is the interaction 
representation of the system coupled to the source. Thus, the physical ob- 
servables of the coupled system will evolve according to the Hamiltonian of 
the isolated system but the states follow the external perturbation, which we 
will refer to as the source. Hence, the expectation value of the observable 
0(x,t) in the exact ground state \G) of H, under the action of the weak 
perturbation H ex is modified as 

(G\0(x,t)\G) -> {G\U-\t)0(x,t)U(t)\G) (3.2) 

where U (t) is the evolution operator in the interaction representation of H?- 
As we have seen, U (t) is given by the time-ordered exponential 

~ f dt'H cxt (t') 
U(t) = Te n J -™ (3.3) 

By expanding U(t) to linear order in the perturbation H cxt , we see that to 
lowest (linear) order the change in the expectation value 5(G\0(x,t)\G) is 

5(G\0(x,t)\G) = \ f dt' (G\[H cxt (t'),0(x,t)]\G) (3.4) 

This change represents the linear response of the system to the external 
perturbation. It is given in terms of the ground state expectation value of 
the commutator of the perturbation and the observable. For this reason, this 
approach is called Linear Response Theory. Notice that in Eq.(3.4), there is 
an ordering of the times t and t': the time t at which the change is observed 
is always later than the time(s) t' during which the external perturbation 
acted, t > t'. Hence, Eq.(3.4) explicitly obeys causality. 

In general, if 0(x,t) is a local observable, H ext (t) represents an external 
source which couples linearly to the observable, 

H cxt (t) = J d 3 x 0(x,t)f(x,t) (3.5) 

To simplify matters, we will assume that the observable 0(x,t) is normal 
ordered with respect to the exact ground state \G), i.e. we will require 
that (G\0(x,t)\G) = 0; in other terms, we are considering operators which 
measure the fluctuations of the observable away from the expectation value. 
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Hence, the linear change (from zero) in the expectation value of the observ- 
able induced by the source f(x,t) is 

(G\0(x,t)\G) f = % -J t dt' J d 3 x' (G\[0(x',t'),0(x,t)}\G) f(x',t') (3.6) 

Thus, the response is linear in the force. The "coefficient" of proportionality 
between the change in the expectation value (G\0(x,t)\G) and the force 
f(x't') defines a generalized susceptibility x(^t^'t') through the definition 

(G\0(x,t)\G) = 'W" = / d 3 x' f dt' x(xt,x't')f(x',t') (3.7) 

J J — oo 

By inspection we see that we can identify the generalized susceptibility 
x(xt;x't') with the retarded propagator, or correlation function, of the ob- 
servable 0(x, t): 

x (xt;x't') ee - % -e(t-t'){G\ [0{x,t),0{x',t')\ \G) (3.8) 

That is, the susceptibility is the retarded Green function of the observable. 
Let f(x,u) be the time Fourier transform of the force f{x,t), 

f(x,t) = jH^ e ™ (3.9) 

The time Fourier transform of Eq.(3.7) is 

(G\0(x,u)\G) = J d 3 x'x(x,x';u)f(x',u) (3.10) 

However, since 
(G\0(x,t)\G) = 

° —e^ [ d 3 x' [ dt' (^f) e^'-'^Gl [0{x,t),0{x',t')] \G)f(x',u) 



2vr J \ h r 

j°°^-e^ J d 3 x'J° dt f (^y^'(G\[0(x,t),0(x',t + t')}\G)f(x',u;) 

(3.11) 
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we find 

(G\0(x,u)\G) = j d s x' (^j j dt'(G\[0(x,t),0(x l ,t + t l )}\G)f(x',uj) 

(3.12) 

Hence, the time-Fourier transform of the susceptibility becomes 

X(oo;x,x') = -- £ dre^(G\[0(x,t),0(x';t + r)]\G) 

= — — I dre^ {G\[O(x,0),O(x' ,r)]\G) (3.13) 

J-oo 

Therefore, we find that the susceptibility x{^^''i u ) is given by 



X(x,x';cu) = -- / dr e lWT (G\[O(x,0),O(x',t)]\G) 



h 



— CO 



(3.14) 



This result is known as the Kubo Formula. 

We just showed that the retarded Green function of the observable 0(x, t), 

V%\x,x') = -i8(t - t'){G\[0{x),0(x')]\G) (3.15) 

is related to the generalized susceptibility, 

1 r°° 

x(co;x,x') = - J dre^V^(x,x') (3.16) 

Hence, 

(G\0(x)\G) = \ J dt' J d 3 x'V r S\x,x')f(x',t') 

= d 4 x'V r S\x,x')f(x') (3.17) 
where x = (x, t) . After a Fourier transform we get 

(G\0(p,u)\G) = ^D^(f,u)f(p,u) (3.18) 
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In Fourier space, the generalized susceptibility xiPi w ) is 

= (GlO(M\G) = 1 ^ Response" 



/I u y A "force" 

, (3-19) 

In conclusion, the response functions of physical interest, and the asso- 
ciated measurable susceptibilities, are given in terms of the retarded Green 
functions of the physical observables. However, as we saw in the last section, 
what we can calculate more directly are the time-ordered Green functions, 
i.e. the propagators of the observables. we will show next that the analytic 
properties of these functions in the complex frequency (or energy) plane 
are such that there is an analytic continuation procedure relating the time- 
ordered and the retarded Green functions. In addition, in practice we really 
are interested in knowing the physical susceptibilities and propagators not 
only in the ground state but also at non-zero temperature. Fortunately the 
analytic continuation will enable us to determine all these functions also at 
non-zero temperature. 



3.2 Finite Temperature 

All physical systems in thermal equilibrium are actually at finite temperature 
T. We will now show that there is a simple and straightforward way to adapt 
the T = methods to account for the effects of thermal fluctuations. We 
will assume throughout that the system is in thermal equilibrium with a heat 
bath at temperature T and at a fixed chemical potential //, i.e. we will treat 
the system in the Grand Canonical Ensemble. The Gibbs density matrix pc 
and the Grand partition Function Zq are 

p G = e-f 3 ( H ~ , Z G = Tr p G = e~^ G (3.20) 

where /3 = l/(kT), H is the Hamiltonian, N is the particle number operator, 
which commutes with the Hamiltonian, [N, H] = 0, Zq is the grand partition 
function and is the grand potential (or thermodynamic potential). 1 
The thermal expectation value, i.e. the expectation value in the the grand 

1 Recall that a chemical potential can only be defined for a conserved quantity, in this 
case the particle number operator N. 
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canonical ensemble, of a physical local observable A(x,t) is given by 

A(1 , t) ) = l^mml ( 3. 21 ) 

Let {| A)} be a complete set of eigenstates of the full Hamiltonian H and of the 
particle number operator N, with eigenvalues {E\} and {N\}. (Notice that 
since we are working in the grand canonical ensemble, necessarily the states 
{|A)} are states in Fock space.) Then, we can write the thermal average as 

A(£ t)\ ~ E,(A|A(f,t)|A) e -^-,N x ) 

Thus, in principle in order to compute a thermal average we will have to 
know first the quantum mechanical expectation value in each state and then 
compute the thermal average. Although this sounds rather laborious there 
are very direct ways to compute the thermal average directly. In what fol- 
lows, as we did before, we will include the chemical potential term in the 
Hamiltonian, H — > H — fiN, which now depends explicitly on the chemical 
potential //. 

To calculate thermal averages it useful to define temperature correlation 
functions. Let A be the Schrodinger operator for an observable. We define a 
temperature variable r, with < r < (3h and the r-dependent operator A{r) 

A(t) = e Hr l h A e~ Hr / h (3.23) 

which is formally obtained by the analytic continuation to imaginary time 
r = it of the Heisenberg operator A(t) 

A(t) = e iHt l h A e - iHt l h (3.24) 

(Notice that in general A\t) ^ A(r)L) 

Likewise, the evolution operator U (t) of a quantum system in the Schrodinger 
picture, now becomes 

U(t) = e - %Ht l h -> U e (r) = e- Hr / H (3.25) 
In particular the Gibbs density matrix Pg(P) is 

p G ((3) = e~P H = U e ((3h) (3.26) 
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Thus a quantum system in equilibrium at temperature T can be regarded 
(formally) as a the imaginary time evolution of a quantum system. In addi- 
tion since the ensemble average is a trace in Fock space, the initial and finite 
states are the same state. In other words, the states obey periodic boundary 
conditions in the imaginary time coordinate r. Thus we will restrict the 
imaginary time variable to the range < r < (3H. 

Much of what we did at T = can be done also at T ^ 0. In partic- 
ular we can similarly define imaginary time ordered products of operators, 
T t (A(ti)B(t 2 )), by analogy with time ordered products. The temperature 
Green function for the electron operator is defined as 



(3.27) 



Gf(x,r;x',T') = -{T^x , r)^(f , r')) / 

Tr (t t (M^t)^1(x',t')) Pg ^ 



Tr pc 

Due to the fact that we are computing a trace it is straightforward to 
show that, as a consequence of the fermionic statistics (and of the anti- 
commutation relations of fermionic operators), the imaginary time propa- 
gator, the temperature Green function, for fermionic operators obeys anti- 
periodic boundary conditions in imaginary time (independently on both r 
and t'). Instead, the temperature Green function for bosonic operators obey 
periodic boundary conditions in imaginary time. Thus 

G T (x, t + f3h; x', t') = TG T {x, t; x', t') (3.28) 

where (— ) holds for fermions and (+) holds for bosons. These periodic and 
anti-periodic boundary conditions have important consequences. In partic- 
ular, since the imaginary time interval is finite, the propagators can be ex- 
panded in Fourier series as 

G T (x - t - t ') = J2 G t(x -x', u n ) e lUn{r - r ') (3.29) 

Periodic and anti-periodic boundary conditions in r and in r' are satisfied if 
the frequencies {u n }, where n G Z, are respectively given by 

w„ = — n, for bosons 

pn 

w„ = — [ n + - J , for fermions (3.30) 

pa \ I ) 
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» ig(n + l/2) 

o ^ Re w 



(a) (b) 

Figure 3.1: Poles on the imaginary frequency axis for (a) fermions and (b) 
bosons. 

These considerations also apply to correlation functions which are always 
bosonic, even in a theory of fermions as they involve operators which are 
typically bilinear in fermions. 

Since the temperature Green function is formally the analytic continua- 
tion of the time-ordered Green function to imaginary time, restricted to the 
imaginary time interval < r < (3h, its imaginary time Fourier transform (or 
series), Gf(x — x', w„) can also be regarded as the analytic continuation from 
real frequency uj to the imaginary frequency axis, restricted to the values iu n 
discussed above for fermions and bosons (See Fig.3.1). 

3.3 Green Functions at T ^ 

We will see now that there is a close relationship between the time-ordered, 
the retarded and the temperature propagators. The same considerations 
apply for the response and correlation functions. The key connection is the 
concept of the spectral function which describes the fluctuation spectrum. In 
particular, for the case of the response functions we will find a connection 
between the fluctuations and the dissipation in the system. 

We will begin with the connection between the fermion propagator, the 
fermion retarded Green function and the fermion temperature propagator. 



Im 00 



2tt 



Re u 



10 



CHAPTER 3. LINEAR RESPONSE THEORY 



To simplify matters I will drop the spin indices (which are easily restored) 
and focus on the time dependence. Except for straightforward sign changes 
the same ideas apply to the boson propagator. Here we will follow closely 
the discussion in Doniach and Sondheimer. 



3.3.1 Spectral Function 

Let us define the fermion correlator G > (a 7 — x',t) 



(3.31) 

Notice that this is a thermal average of a product (not necessarily ordered!) 
of Heisenberg operators. 

Let {|A)} be a complete set of eigenstates of the Hamiltonian in Fock 
space. Completeness means that the identity operator in Fock space can be 
expanded as 

1 = $>>(A| (3.32) 

A 

Then, we can write the correlator as 



^(x,t) ^(x',0)) = -^-Y (AI^O?' )e-P H \\'){\'\e iHt/n iP(x)e- iHt/R \X) 

= _L£V^v(A|^t(£' )\X')(\'\,l;(x)\X)e i ^'- E ^ h 

(3.33) 



U XX' 



where 

A 

is the (grand) partition function. 

Let Ji(x — x';u) be the time Fourier transform of {ip(x, t) ift{x", 0) ), 



Ji(f-x';w) = / U(x,t)^(x\0)\e iuJt dt 

J —CO 

^^e-^(A|^(f')|A')(AV(f)|A) 2vr5(^-^ + c) 

G XX' 

(3.35) 
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J\{x — x'; oj) is the spectral function of the time correlator (ip(x, t) ip^(x', 0)). 

Similarly we can also write the correlator G < (x—x / , t) = (^p^(x', 0)ip(x, t) 
and its time Fourier transform J 2 (x,x';uu): 

J 2 (x-x';u) = f U\x',0)^(x,t)\e iujt dt 

J —CO 

= j- E e ^ Ex W f ( f 1 y ) ( A >( f ) I A > 2?r 5( Ex '~ Ex + co) 

G XX' 

(3.36) 

Using the properties of the Dirac 5-function we get 

J 2 {x-x'-u) = —Ye-^ E * + ^)(\\^{x')\\'){\'\i}{x)\\) 2n5( Ey ~ Ex + u) 
= e'^J^x-x'-u) (3.37) 



Notice that at T = 0, only the state |A') = \G) (the ground state) survives 
in sum over A'. Thus, at T = 0, J±(x — x';u) vanishes for to < 0, while 
J 2 (x — x'; to) vanishes (at T = 0) for u > 0. 

3.3.2 The Retarded Green Function 

Let us consider now the retarded fermion Green function G R (x — x',t) (at 
finite temperature), 

G R (x-x',t) = -ie(t)^{ip(x,t),^(x',0)}^ 

= -i9(t) {{^(x,t) ^(x',0)) + (^(£',0)^,0 )) 
= -i6(t) (g>{x -x \t) +G K (x -x ',*)) 

(3.38) 

where we have used anti- commutators because we are dealing with fermion 
operators. For a bosonic operator we must replace the anticommutator by a 
commutator. 

For a translationally invariant system we can clearly express the retarded 
Green function in terms of the spectral functions J\(k,u) and J 2 (k,u), the 
space Fourier transforms of the functions Ji(x,x'; oj) and J 2 (x,x'; oj) (where 
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k is the wave vector) defined above. To do that we will note that for the case 
of translationally invariant systems the Hamiltonian also commutes with the 
total linear momentum operator P, [P, H] = 0, and hence the states {|A)} 
can also be chosen to be eigenstates of the linear momentum P, with eigen- 
values {-P\}- If we now recall that since P is the generator of infinitesimal 
translations in space, we can also write 

^(f) = e ip -* /h i>{<S)e- tp - i,h (3.39) 

— * 

Hence, in Fourier space, the spectral function Ji(k,u) becomes 

G XX' 

(3.40) 

where P\ and Py are the total linear momentum of the states |A) and |A'). 
Notice that Ji(k,u) is a real function. 

By the same argument used above we also get 

J 2 (k,u) = e- f " iu J 1 (k,u) (3.41) 

Then, the retarded fermion Green function at T > is 

G R (x,t;x',0) = -i9(t) ^ (jx&u/RJa&o/)) e^*"*')-") 

J (2vr) >J 27r V J 

(3.42) 

Its Fourier transform G R (k,u) is 

G R (k, u) = — ( 1 + e-^') JliKuJ,) (3.43) 
y.oo 2tt V J u-u' + iri 

(where i] — > + ). This function has poles in the lower half plane Imcu < 0. 
We will now use the fact that J\(k, u) is a real function and the property 

lim — !— — V- - m5(x) (3.44) 

»7^0+ X + If] X 
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where x is real and 

V- = lim — ^— (3.45) 
is the principal value of 1/x, to show that 

/oo j / . . 
(l + e _/3W ) Ji(jfe, w') <J(w - w') 

= ~(1 + e-^)Mlz, U ) 

(3.46) 

Hence we obtain 

Wlv) = - + (3.47) 

For a bosonic operator, and in particular this will apply to the general- 
ized susceptibilities, the corresponding relation between the spectral function 
Ji(k,uj) and the imaginary part of the Green function is 



3.3.3 The Time-Ordered (Feynman) Fermion Propa- 
gator 

Let us consider now the Feynman fermion propagator (also at T > 0) 

G F (x -x',t) = -i(T(ip(x, t)ip\x', 0)\ (3.49) 



where we use a fermion time-ordered product. It can also be expressed in 
terms of the spectral function Ji(x — x';oj) and of its Fourier transform (in 
space) Ji(k,u): 



/OO J / 
—j x {x-x\jy- iwlt 
-oo 2vr 



oo 

00 duo' 



+i0(-t) I ^J 2 (x-x',co')e- tu ' t 

l-oc ^ 



(3.50) 
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and its time (and space) Fourier transform 



OC 



27r 1 w — u/ + «?7 oj — u' — it] 



00 doj' T . r ,. f 1 e"^' 
Ji(fc,a; ) 



■DO 



27r ' [ w — + irj oj — oj' — irj 



(3.51) 



which has poles both in the upper and in the lower half of the complex oj 
plane. Thus, contrary to the retarded Green function, the Feynman propa- 
gator is not analytic on either half plane. 
Since Ji(k,u>) is real, we get 

Re G R (k, oj) = Re G F (k, oj) = V f° ( 1 + e'^') ^tl^l d ^L (3.53) 

J -00 v ) oj -oj' In 

and 

lmG R (k,oj) = ~{l + e - fihu ) Ji(k,u) 

1m G F {k,u) = -~{l-e- phw ) Ji(k,oj) (3.53) 

which imply that the real and imaginary parts of the retarded and time 
ordered functions are related by 

ReG R (k,oj) = ReG F (k,oj) 

lmG R (k,oj) = coth(^)lmG F (k,oj) (3.54) 

At T = 0, 

Im G R (k, oj) = sign(w) Im G F {k, oj) (3.55) 

Finally the real and imaginary parts of the retarded (and time-ordered) Green 
function are related by the Kramers- Kronig relation 

t. ^n,7 \ ^ f°° du'lmG R (k,oj') 

ReG R (k,oj) = -V — — '- 3.56 

J_n 7T OJ-OJ' 
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3.3.4 The temperature Green Function 

The temperature Green function Gt(x—x', t) is the time-ordered propagator 
in imaginary time r, 

G T (x - x', r) = -(T T (iP(x, t)^(x', 0))^ (3.57) 

where T T is the imaginary time ordering operator, and < r < f3H. Hence 

G t (x-x',t) = -9(r){ij(x,T)^(x',0)) + 9(-T)^(x',0)ij(x,T)) (3.58) 

Once again we expand the thermal expectation values in terms of complete 
sets of states {|A)}, to obtain 



AA' 

(3.59) 

Since the temperature (or imaginary time) fermion Green function is anti- 
periodic in time, Gt(x — x',t) = —Gt(x — x',t + (3h), it can be expanded 
in Fourier series with coefficients 

rfih 



iP(x,r)^(x',0)) = -^V e-^'(X\^\x')\X'}{X'\^(x)\X}e T(E ^~ E ^ /h 
1 Zg 7t 

i>\x',0)4>(x,T)) = ^- Ve-^(A|^ t (f / )|A / )(A / |^(f)|A)e T ^'-^^ 
/ Z G ^- 



i r h 

G T (x-x',uj n ) = — dre-^ T G t (x-x',t) 
P* 1 Jo 

1 



1 r pn 

- / dTe- tu " T (^(x,r)^x')) 
Jo 

Y,W\x')\X')(X'\^{x)\X) 1 " + < V 



P hZG AV {Ey -E x - ihWr, 



(3.60) 



where u; n = g(n + 1/2). 
In Fourier space we find 



_ (7 , 1 f°° du' / J 1 (k,tU>) + J 2 (k,tu') 



(3h J_ 00 2n \ —iuJ n ~ oj' 

f3H J.^ 2-K \ J -iu n - U)' 



(3.61) 
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Im to 

i||(n + l/2) 



' pn 



Re oj 



Figure 3.2: Analytic continuation to real frequency at finite temperature. 



Thus, the temperature Green function Gr(k, oo n ) can also be determined from 
the spectral function Ji(k,uu). 

Furthermore, if we now compare this result with the analogous expression 
for the retarded Green function G R (k,u) we see that 



G R (k,-ico n ) = -PhG T (k,uj n ) 



(3.62) 



— * — * 

Hence, Gn(k,ui) is the analytic continuation of —(3hGT(k,u> n ) to a patch of 
the complex plane, that includes the imaginary axis: 



—iuin^-uj + ir] =>• —/3hGT(k,uj n )^-G R (k,uj) 



(3.63) 



Since GV(&, w n ) is known only for a discrete set of frequencies on the imag- 
inary axis, {ioo n }, this procedure amounts to an analytic continuation from 
the imaginary frequency axis to the whole complex frequency plane, and thus 
to the real axis. (See Fig. 3. 2) 

Thus we have now a direct way to compute these functions 

1. We first compute the temperature Green function GT(k,u n ) 

2. Next, by analytic continuation we find the retarded Green function 
G R (k,u) 

3. The spectral function Ji(k,u) is obtained from the imaginary part of 
the retarded Green function. 
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4. The time-ordered and the retarded Green functions are both deter- 
mined in terms on the spectral function 



3.4 Dissipation and Response 

We will now extend these ideas to response functions. We will find along the 
way a very important result known as the fluctuation- dissipation theorem. 

We will begin by considering first a very simple problem. Imagine that 
we have a classical linear harmonic oscillator, represented by a coordinate 
x(t). The oscillator has mass m and natural frequency uq. The oscillator 
is in thermal equilibrium with a bath represented by a set of random time- 
dependent "internal forces" F int (t), which represent the collisions between 
the degrees of freedom of the bath with the oscillator. Thus in equilibrium 
the classical equation of motion of the oscillator is 

+ mcu 2 x(t) = F int (t) (3.64) 

We will now imagine that we act with an external force F ext (t) on the os- 
cillator. 2 The average non- equilibrium displacement (averaged with respect 
to the random internal forces) which we will denote by (x(t)) n . e ., is linearly 
related to the external force by an expression of the form 

POO 

(x(t)) n . c .= x(t,t')F cxt (t')dt' (3.65) 



The equation of motion in the presence of the external force is 

d 2 x(t) 

m ~d^ + mw oM*) = *int(f) + ^cxt(t) (3.66) 
so that on average we have 
/ d 2 x(t) \ 

m \—^T-) nc + mu 2 {x(t)) n . e . = (F int (0) n .e. + F cxt (t) (3.67) 

where we have used the fact that, in the presence of the external force, the 
average effect of the internal forces, which we define as friction, does not 



2 Notice that in the presence of the external force, the system (the oscillator) is no 
longer in equilibrium. 
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vanish, (Fi nt (t)) n , e , 7^ 0. In simple phenomeno logical models we normally use 
the "constitutive relation" which asserts that the friction force is a linear 
function of the velocity of the oscillator 

(F int (t)) Q . c . = - mi (^±}^ (3.68) 
where 7 is the friction constant. The equation of motion now is 

m (^^)„c +m7 (^)ne +mW o(^)>n.c. = F cxt (t) (3.69) 

Hence, x(M') satisfies 



m- 



d 2 x(t,t>) , dx(t,t') 



+ m l \l + m "oX(t, = S(t - t') (3.70) 



dt 2 ' ' dt 



which is to say that x(t, t') = x(t — t') is the retarded Green function for this 
non-equilibrium system. 

It is straightforward to show that x(t — t') has the integral representation 

X(t-t') = / = : — 3.71 

where the integration runs along the real axis on any path in the upper half 
plane. Thus, the real and imaginary parts of its Fourier transform, x'( w ) an d 
x"(w), in this simple problem are, respectively, given by 

x-'m = r " g :r 2 — * (3^72) 




X» = -r (3-73) 



(u, 2 -u, 2 ) 2 + ( 7 c) 2 ] 



As usual, x'( w ) an d x"( w ) obey the Kramers-Kronig relation. 

Using Cauchy's Theorem we can define a function x( z ) where z is a 
complex number, 



dwxH rdw X "(co) 



/- (3.74) 
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where z is in the upper half plane, Imz > 0. By explicit calculation we find 

X(z) = -- ^— (3.75) 

m z 2 — Uq — i^z 

which does not have a pole on either half plane but instead a branch cut 
along the real axis, with discontinuity 2%"(u;). Notice that the branch cut 
arises only in the presence of a finite damping coefficient 7. In other words, 
it is due to the existence of friction forces. 

We will now show that x"i u ) represents the dissipation. Let us compute 
the work done by the external force, per unit time, 



dW /dx{t)\ 
— if 



dt ext \ dt I n 

duj' 



FUt) f dt'^(t-t')F cxt (t') 



/OO J / poo 
^ e-*-* F cxt (t) {-iJ) X {u') / e-'*'F ext (^) 
-00 ^ J -00 

(3.76) 

For a "monochromatic" force F ext (t) = Re (J^xte^*), the power dissi- 
pated, averaged over one cycle, is 



dW cu 



(H ^/H^W (3-77) 

where we used that x"( w ) = ~x"(~ uj )j which holds for this simple model 
but it is more generally true. Hence, we see that the dissipated power is 
proportional to x"( w )j ^- e - to the discontinuity across the branch cut. In 
particular this simple model suggests that there is a direct connection be- 
tween dissipation and susceptibilities. This result is actually quite general 
and we will identify the imaginary part of the response functions with the 
dissipative response. 



3.5 The Fluctuation-Dissipation Theorem 

Let us now return to the discussion of the response functions. As we saw 
above, even in Fermi systems the response functions are related to correlators 
of bosonic operators (e.g. bilinears of Fermi operators). In particular, for a 
general bosonic operator 0(x, t), the generalized susceptibility x(%, t; x', t') is 
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simply related to the retarded Green function V ret (x, t; x', t') for the operator 
0(x,t): 

X (x,t-,x',t') = ^ ct (x,t-x',t') (3.78) 

This relation is also correct at finite temperature T. Thus, just as in the case 
of the fermion correlators, we will need the connection between the retarded, 
time-ordered and temperature correlation functions. These relations are very 
similar expect for minor changes due to Bose statistics. We will not repeat 
all the arguments and we will instead just write down the important results. 
Just as in the fermionic case, the key concept is the spectral function. 

Thus, for a general bosonic operator 0(x,t) we define the correlators 
£>>(£-£' ,t), £><(£-£ ',t), V rct (x-x',t), V F (x-x',t) and V t (x-x',t), 
by the thermal averages 



Z>> (£-£',*) = (O(x,t)O(x',0) 
D<(£-£',t) = (o(x' ,0) 0(£,t)) 



V vct (x -x ',*) = -i0(f)^[O(£',f),O(£,O)] 
V F (x-x',t) = -i(r(O(x',t),O(x,0)) 
V t (x-x',t) = (t t (O(x',t),O(x,0)) 



(3.79) 

— * — * 

The spectral functions J\{k,u) and J2(k,u) are once again the space and 
time Fourier transforms of the correlators D > (x — x',t) and T> < (x — 
respectively: 

— * — * 

Ji(k,u) = V > (k,co) 

J 2 (k,u) = P<(jb,o;) (3.80) 

which obey 

J 2 (£,w) = e"^ Ji(k,u) (3.81) 

However, since the operators 0(x,t) are bosonic (that is, they obey com- 
mutation relations), the temperature Green function obeys periodic bound- 
ary conditions in imaginary time, T>t(x — x',t + f3h) = T>t(x — x',t). This 
condition, which implies that the Matusbara frequencies now take the val- 
ues u n = ||n, where n G Z, has important consequences. In particular, 
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paraphrasing what we did for the Fermi case, we now find that the spec- 
tral function J\(k, oo) and the imaginary part of the Fourier transform of the 
retarded Green function Im V iet (k, u>) are related by 

Ml") = - l _l_ ph J mV^{kM (3.82) 

In particular, this relation implies that the dissipative component of the 
dynamical susceptibility x"(k,u) is proportional to the spectral function 
Ji(k,u), 

x >\k,co) = ~(l-e-^).J 1 (k,cu) (3.83) 

This relation is known as the Fluctuation-Dissipation Theorem. 

On the other hand the retarded and time ordered Green functions are 
now related by 

ReD rct (k,oo) = ReV F (k,u) (3.84) 
lmD TCt (k,co) = tanh(^)Im£> F (£,w) (3.85) 

Also, just as in the fermionic case, the real and imaginary parts of the re- 
tarded correlation function are related by the usual Kramers- Kronig relation 

B»V~& U ) = -V f^jSZ^O (3.86) 



— oo 



Finally, the temperature Green function V T (k,cu n ), where u n = is 
related to the spectral function by the integral transform 



IMj U) = 4f£(i- e -~)4<H (3. 



87) 



Thus, here too, all the correlation functions of interest can be found by the 
same analytic continuation procedure we already discussed for fermions. The 
bosonic nature of the correlation functions follows form (a) the physically 
correct spectral functions and (b) the temperature-dependent prefactors in 
the integrand of this equation. 
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3.6 Perturbation Theory for T > 

We will now discuss briefly perturbation theory at finite temperature T. The 
approach is very similar to what we did at T = 0. 

Let is consider a system with Hamiltonian H = Ho + Hi. We will be 
interested in the fermion temperature Green function Gt(x — x', r) 

G T (x -x',t-t') = (t t (^(x, t)^{x', t')) (3.88) 

for a system described by the Hamiltonian H. 

We now introduce a representation analogous to the interaction represen- 
tation but for T > 0. Let U (r, r') be the operator 

U(T, T') = e rHo/h e -(r-r')(H +H 1 )/h e -r'H /h ^gg) 

Notice that although this operator is not unitary, it satisfies the group prop- 
erty 

U(t,t')U{t',t") = U{t,t") (3.90) 
and it obeys the (imaginary time) equation of motion 

-H^-(t,t') = Hi(t)U(t,t') (3.91) 

where 

H^t) = e rHo/n Hie -rH /h ^ 92 ) 
In particular the full density matrix is given by 

e^ 11 = e~ l3Ho U(/3h,0) (3.93) 

It is straightforward to show that the temperature Green function of the 
full system is given by the following expression in terms of thermal averages 
of the unperturbed system: 



T T (U((3h,0)^(x,T)^(x',r')) 

G t (x-x',t-t j ) = + -. r ^ (3.94) 

U(f3h,0) 





Here we used the notation 



Tr (e~^ H °A) 
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where A is an arbitrary operator. 

Notice that the operator U(/3h,0)) plays here a role analogous to the 
S'-matrix. Hence, using the same line of reasoning that we used for the 
S'-matrix, we can write 



Thus we can carry out perturbation theory much in the same way as we did 
at T = 0. The only missing ingredient is a generalization of Wick's Theorem 
for T > 0. It turns out that this works exactly in the same way provided that 
the contraction of an arbitrary pair of operators A{ji) and B(t 2 ) is identified 
with imaginary time ordered thermal expectation values of the operators with 
respect to H : 



which is a thermal contraction. Notice that we are using the same notation 
for the thermal contraction that we used at T = although it now has a 
different meaning. 

Wick's Theorem simply states that the thermal average of the product of 
any number of operators is equal to the sum of the products of all possible 
pair-wise thermal contractions. Consequently, the structure of the perturba- 
tion series is the same as the one we found at T = 0. In particular we will 
also find that disconnected diagrams, which at T = are called the "vac- 
uum diagrams", also cancel out exactly at T > 0. They only survive in the 
computation of the grand partition function. Thus, all of the Feynman rules 
used at T = carry over to T > with the only change that the factor of 
—i/h now simply becomes —1/h. All other rules remain intact, except that 
the propagators are now temperature Green functions, which obey specific 
periodic and anti-periodic boundary conditions in imaginary time. 



3.7 The Electrical Conductivity of a Metal 



We will now consider the response of an electron gas to weak external elec- 
tromagnetic fields Ap(x). The formalism can be generalized easily to other 





(3.96) 




(3.97) 
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systems and responses. In particular, we will consider the electrical conduc- 
tivity of a metal. 

There are three effects (and couplings) that we need to take into consider- 
ation: a) electrostatic, b) diamagnetic (or orbital) and c) paramagnetic. The 
electrostatic coupling is simply the coupling to an external potential with 
H cxt given by 

#cxt = Yl / ^ x e( t>( x ^) ij>l(x)Mx) = J d * x Jo(x)Mx) (3-98) 
c=t,4- 

where = A is the scalar potential ( or time component of the vector po- 
tential Ap). The diamagnetic coupling (or orbital) follows from the minimal 
coupling to the external vector potential A. The kinetic energy term -f^in 
is modified following the minimal coupling prescription to become 



H kin (A) = I d 3 x^- f v + ^0*)) ^( x ) ■ (v - ^/{x) ) ^(x) (3.99) 



which can be written as a sum of two terms 



H kin (A) = H kin (0) + H cxt (A) 



(3.100) 



where i^kin(O) is the kinetic energy term of the Hamiltonian in the absence 
of the field, and H ext (A) is the total perturbation, i.e., 



H ext (A) 



d 3 x 
d 3 x 



J (x)A (x) - J(x) ■ A(x) - ^—^A 2 (x)J (x) 



J d 3 x [j^x)A^(x) - ^-^A 2 (x) Jo(rr)] (3.101) 



Here Jo(x) = ep(x) = e^o-Vo-Vv is the local charge density, and J(x) is the 
gauge-invariant charge current defined as 



J(x) 



ieh 
2mc 

ieh 
2mc 



Y,WAx)D^{x) - (d^{x)) ] 4, a {x)] 



(3.102) 
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— * — * — * 

where D = V — ij^A is the (spatial) covariant derivative. We are using here 
a 4-vector notation J M = (J , J) and J^A^ = JqAq — J ■ A. 

Clearly J(x) is the sum of the two terms, one which represents the mass 
current and the diamagnetic term, -^Atp^ip. We can write the total pertur- 
bation, including the scalar potential A , if we write 

H cxt = J d 3 x[J,(x)A^x) - ^— 2 A\h{x)\ (3.103) 

Finally, we can also consider a paramagnetic coupling to the spin degrees 
of freedom which has the Zeeman form 

^Zeeman = J fa . ^S^(x) (3.104) 

where g is typically of the order of the Bohr magneton hb and S = | for spin 
one-half systems. We will discuss the Zeeman term later in the context of the 
magnetically ordered states. For now will focus on the purely electromagnetic 
response of a general interacting Fermi system. 

A straightforward application of the Linear Response formulas derived 
above yields an expression for the current (J^)' in the presence of the per- 
turbation. 

«*)>' = (W)G- l j-f_ dt' Jd 3 x'(G\[J„(x'),J tM (x)]\G)A„(x') + ... 

(3.105) 

This formula suggests that we should define the retarded current correlation 
function T> r ^(x, x') 

V^(x,x') = -iQ{x Q -x' G )(G\[J,(x),J u {x')]\G) (3.106) 

The induced current (Jk)md 

(J fe )ind = (JkY - (jk) G (3.107) 

(where j k is the "paramagnetic" component of the current) has a very simple 
form in terms of V r ^;(x, x'), namely 

= \J d 4 x f ±{vX(x : x>)Mx')-V*(x,x')A (x')} 

e —J Q {x)A k {x) + ... (3.108) 



J k (x) 

ind 



mc* 
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3.7.1 The dielectric tensor and the conductivity tensor 

The induced current has two important properties: a) it is conserved, namely 
the the induced current and the induced density satisfy a continuity equation 

^f + V-4d = (3.109) 

and b) it is gauge invariant. This means that we should be able express the 
induced current directly in terms of gauge invariant quantities such as the 
electric and magnetic fields. 

Since (J^x))^ is gauge invariant, we can compute its form in any gauge. 
In the gauge A = the spacial components of ( J M (x))i n d are 

(Jk(x))- md = -^A k (x) + - / d 4 x'Vlf(x - x')A e (x') + ... (3.110) 
mc z n J 

— * — » 

In this gauge, the external electric field E ext and magnetic field H are 

E cxt = -d A H = VxA (3.111) 
Now, in Fourier space, we can write 

e 2 p 1 
( J k (p, w)) ind = -^A k (p, u) + -^D r k f(p, u)A e (p, u) 

1 p 2 n \ F! cxt 

VlUp,u)-^-5 k A^{p,uj) 



fi mc z J iu 

(3.112) 

This expression is almost the conductivity. It is not quite that it since the 

— * — * — * 

conductivity is a relation between the total current J = J m( \ + J ext and the 
total electric field E, instead of one between the induced current and the 
external electric field. In order to take these electromagnetic effects into 
account, we must use Maxwell's equations in a medium which involve the 
fields E, D, B and H 



v .D = p vx£ = - 9J! 
V-5 = vx^f + J 



dt (3.113) 



where 

B = H + M E = E cxt + E md (3.114) 
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Here M and E md are the magnetic and electric polarization vectors. In 
particular 

J md = d t E md (3.115) 

and 

d t D = d t E + J md (3.116) 
Linear Response theory is the statement that D must be proportional to E 

D j = e jk E k (3.117) 

where Sjk is the dielectric tensor. Since E and E ext satisfy similar equations 

-V x y x E = d 2 E + 8 t J 
-V x v x E cxt = d 2 E ext + d t J ext 

(3.118) 

and v x V x -E = v(V ' -^0 ~ V 2 -^? we can write, for the Fourier transforms, 
the equations 

PiPjEj(p, co) - p^Ei(p, co) = -u 2 Ej(p, u) - iuJi{p, uj) 
p iPj Ef (p, co) - tfEf* (p» = -co 2 E^ (p» - iuJ?* (p, u) 

(3.119) 

Thus, we get 

PiPjEj(p, co) - fEi(p, co) + uJ 2 Ei(p, co) = -iuJl nd (p, co) 

+ p tPj E] xt (p, co) - i?Er\p, co) + co 2 Er\p, co) 

(3.120) 

Since we showed above that 



-icoJ^ip, co) = ( 5^ - ±Df{p, co) ) Ef{p, co) (3.121) 



we conclude that 

(PiPj -ffiij + co 2 5ij) Ej(p,co) = 

(Sij^ - lv*(p, co) + p iPj - p% + co% 3 ) E] x \p, co) 

(3.122) 
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In matrix form, these equations have the simpler form 

{p®p-fl + u 2 I)E = (^1 - \v R +p®p-fl+ cu 2 l] E cxt (3.123) 

\mc z n J 

This equation allows us to write -B ex t i n terms of E. 

Hence we find that the induced current is related to the total field by 

\ mc z J mc z a 

(p®p-j?I + oo 2 l) E 

(3.124) 

It follows that the conductivity tensor a is 

iuaftu) = (\v^{p^)-^l\ + (\v^{p,u)-^l) 
\n mc z J \n mc z J 

_ lv IC \p, oo)+p®p-fl + u 2 I]- 1 
mc z a 

\ n mc z J 

(3.125) 

— * — * 

Also, since D = eE, the dielectric tensor e is 

e = I + -a (3.126) 

to 

We conclude that both the conductivity tensor and the dielectric tensor can 
be determined from the retarded current and density correlation functions. 

3.7.2 Correlation Functions and Conservation Laws 

In the problem discussed in the previous section, we saw that we had to 
consider a correlation function of currents. Since the currents are conserved, 
i.e., <9 M J M = 0, we expect that the correlation function V^x, x') should obey 
a similar equation. Let us compute the divergence of the retarded correlation 
function, d^V^ t (x,x'), 

dX t (x,x') = d*[-iQ(x -x' )(G\[J»(x),r(x')]\G) (3.127) 
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Except for the contribution coming from the step function, we see that we 
can operate with the derivative inside the expectation value to get 



apzM = -i(dze(x -x' )){G\[nx),r(x')]\G) 
-ie(x -x' )(G\[d*j»(x),r(x>)]\G) 

(3.128) 

The second term vanishes since J^{x) is a conserved current and the first 
term is non zero only if fi = 0. Hence we find 

dpZ(x,J) = -i5{x,-x' Q )(G\[J {x),r{x>)]\G) (3.129) 

which is the v.e.v. of an equal-time commutator. These commutators are 
given by 



(G\[J (x,x ),J (x',x )]\G) = 



• 2 

ie 



(G\[J°(x,x Q ),J i (x',x Q )]G) = 



mc 2 



(3.130) 



Hence, the divergence of V r ^ is 

d£v$M = ^-M^ ~ ; = ( 3 - 131 ) 

and 

d:,Vl*{x,x>) = -^dl[8\x-x')(p{x'))\ ■ d^(x,x') = (3.132) 

Notice that the time-ordered functions also satisfy these identities. These 
identities can be used to prove that (J md ) is indeed gauge- invariant and 
conserved. Furthermore, in momentum and frequency space, the identities 
become 



-iu>V$&u>)-i P kD%(p,u>) = 
-iuVlt(p,u)-i Pe Vlf(p,uj) = - 
-iu>V$(p,u)-ip k l>$(p,u>) = 

2 

-iuV k e *(p,uj)-ip e V k f{p,uj) = ~^ip k 

TYIC 



mc z 



(3.133) 



30 



CHAPTER 3. LINEAR RESPONSE THEORY 



We can combine these identities to get 

2 

uPVftlp, u) - PtPkVgiP, w) = -^4p 2 (3.134) 



mcr 



Hence, the density-density and the current-current correlation functions are 
not independent. A number of interesting identities follow from this equation. 
In particular if we take the static limit co — >■ at fixed momentum p, we get 



,2 



Hmp^(p, co) = ^Lp 2 (3.135) 



provided that ]im u ^ V(^(p,u}) is not singular for p ^ 0. Also from the 
equal-time commutator 



we get 



ip 2 

(G|[J fc (f,x ), J (f ,x )]|G) = —d x k (6(x -af){p(x))) (3.136) 

mr 

Jim dZVlf(x,x>) = — 2 vl (6(x-x')(p(x))) (3.137) 
x' — >xo mc 

If the system is uniform, (p(x)) = p, we can Fourier transform this identity 
to get 

The conservation laws yield the alternative expression 

= (3^39) 

This identity is known as the f-sum rule. 

If the system is isotropic, these relations can be used to yield a simpler 
form for the conductivity tensor. Indeed, if the system is isotropic, T>lf(p, co) 
is a sum of a longitudinal part V™* and a transverse part V 1 ^ 

V%(p, co) = Vf(p, co)^ + Vf(p, co) - 5^ (3.140) 

Thus, we get a relation between and the longitudinal part £>j| Ct 

co 2 V^(p, co) - p 2 Vf(p, co) = -^p 2 (3.141) 
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Hence 



and 



oj 2 \ 11 mc z 



lim Vf(p,u 



e 2 p 



(3.142) 



(3.143) 



for all p. 

The conductivity tensor can similarly be split into longitudinal cry and 
transverse a± components 



"13 



PiPj . 



P 2 ^' 



(3.144) 



We find 



and 



(Til = 



VT - 



e 2 p 
mc 2 



ret 



i + 



- + w 2 - p : 



(3.145) 



(3.146) 



These relations tell us that the real part of cry is determined by the imaginary 
part of P|| Ct . Thus, the resistive part of cry (which is responsible for dissipation 
in the system) is determined by the imaginary part of a response function. 
This is generally the case as it follows from the fluctuation-dissipation theo- 
rem. 



3.7.3 Consequences of the /-sum rule 

We will now assume that there is a stable collective mode, i.e., a plasmon 
branch. Therefore the retarded density- density correlation function T>q must 
have a pole at the plasmon dispersion, T\w v \{p): 

V ™ M = J=^) + ^ (3 - 147) 

where A(p, uj) is the residue, and [Dqq\ reg is regular (analytic) near the plasma 
branch. In order to enforce the condition that T>q Q is the retarded correlation 
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function we will perform the analytic continuation u — > cu + ie. This analytic 
structure of the correlation function holds for oo ps |oy(p)|. 

We will now see that the sum rule determines the residue A(p,u). How- 
ever, since [^ > oo] re is non-singular, we have 



/+oo 
duu [V§ ) ] mg (u,p) = (3.148) 
-oo 



since we can close the contour on the upper half plane without enclosing any 
singularities. Thus, in order to satisfy the /-sum rule we must demand that 

+ °° . A{p,u>) e 2 p 2 

du ioj- — = — - p (3.149) 



(to + is) 2 - oj 2 e (p) mc 2 



We will now close the contour on the pole and compute the residue. The 
pole is at w = ±|a> p ^(p)| + ie, and the residue is 

residue: = — nA(p, ±|oy(p)|) (3.150) 

which leads to the a condition for the amplitude A(p,u>): 

2 

-7T(A(p,+\uj pe (p)\) + A(p,-\u pe (p)\)) = ^p 2 (3.151) 

mcr 

Hence, 

A = ~ ^ P 2 (3.152) 

(typically A(u) = A*(-u)). 

Thus we conclude that for frequencies close to the plasma branch the 
retarded density-density correlation function must have the form 

Notice that the residue is determined by the exact density and by the bare 
mass. This is a consequence of Galilean Invariance. 



